#include "NewtonInterpolation.hpp"

using namespace std;
std::vector<double> x = {0, 1, 2, 3};
std::vector<double> y = {1, M_E, M_E*M_E, M_E*M_E*M_E}; 
NewtonInterpolation F1(x,y);
double x_value = 1.5;

double f  (double x) { return pow(M_E,x) ;};
NewtonInterpolation F2(x,f);

int main() {
cout << F1(x_value) << endl;
cout << F2(x_value) << endl;
}
